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The computation of radiative energy loss in a finite size QCD medium with dynamical constituents 
is a key ingredient for obtaining reliable predictions for jet quenching in ultra-relativistic heavy ion 
collisions. We here present a theoretical formalism for the calculation of the first order in opacity 
. . . radiative energy loss of a quark jet traveling through a finite size dynamical QCD medium. We show 

On ' that, while each individual contribution to the energy loss is infrared divergent, the divergence is 

' naturally regulated once all diagrams are taken into account. Finite size effects are shown to induce 

, a non-linear path length dependence of the energy loss, recovering both the incoherent Gunion- 

Bertsch limit, as well as destructive Landau-Pomeanchuk-Migdal limit. Finally, our results suggest 
{^j^' a remarkably simple general mapping between energy loss expressions for static and dynamical QCD 

^ media. 
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I. INTRODUCTION 

^ , . . 

^ ' Suppression pattern of high transverse momentum hadrons is a powerful tool to map out the density of a QCD 

^ . plasma created in ultra-relativistic heavy ion collisions [1-3]. Since suppression (called jet quenching) results from the 
energy loss of high energy partons moving through the plasma [4-7] , reliable theoretical predictions for suppression 
, require reliable energy loss calculations. 
^ ' The medium-induced radiative energy loss is, in most studies, computed by assuming that the QCD medium consists 
of randomly distributed static scattering centers ("static QCD medium"). However, in reality, constituents of the 
■ medium are dynamical, and we recently showed that inclusion of dynamical QCD medium effects are important in the 
radiative energy loss calculations [8, 9]. Furthermore, calculation of the energy loss has to be performed in a finite size 
QCD medium, since the size of the QCD medium created in both RHIC and LHC is finite. While methods for energy 
' loss calculation have been developed for infinite optically thick dynamical QCD medium, no such approach exists for 
, finite optically thin medium. However, it is of crucial importance to develop the energy loss formalism for the case of 
finite size optically thin medium, in order to make reliable predictions applicable for the range of parameters relevant 
for RHIC and LHC experiments. 

This paper develops theoretical formalism for the calculation of the radiative energy loss in a finite size dynamical 
QCD medium, while main numerical results of the model are presented elsewhere [9]. Our work is based a novel 
approach, where two hard thermal loops arc implemented within a finite size QCD model initially introduced by 
' Zakharov [10]. The computations are presented in the appendices, and include analytical calculations of 24 Feynman 
diagrams, each of which is individually infrared divergent. However, this divergence is naturally regulated when 
all the diagrams are taken into account. Furthermore, we obtain an explicit analytical expression for the energy 
loss in dynamical medium, which can be directly compared with the equivalent expression obtained in static QCD 
medium. Finally, we discuss the obtained analytical result in the context of several qualitative features: i) Recovery 
of the static approximation for asymptotically high values of energy, ii) Transition of the thickness dependence from 
Gunion-Bertsch (GB) [11] to Landau-Pomeanchuk-Migdal (LPM) [12] limit, and Hi) Possible extension of the results 
to higher orders in opacity. 



II. RADIATIVE ENERGY LOSS IN A DYNAMICAL QCD MEDIUM 

In this Section we outline the computation of the medium induced radiative energy loss for a heavy quark to first 
order in opacity. We consider finite QCD medimir of size L and assume that the heavy quark is produced inside the 
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medium at time xq = 0, at the left edge of the medium, travehng right. 

Medium induced radiative energy loss is caused by the radiation of one or more gluons induced by coUisional 
interactions between the quark of interest and partons in the medium. The energy loss rate can be expanded in 
the number of scattering events suffered by the heavy quark, which is equivalent to an expansion in powers of the 
opacity. For a finite medium, the opacity is given by the product of the density of the medium with the transport 
cross section, integrated along the path of the heavy quark. The lowest (first) order contribution corresponds to one 
collisional interaction with the medium, accompanied by emission of a single giuon. 

We compute the medium induced radiative energy loss for a quark jet to the first (lowest) order in number of 
scattering centers. The finite size medium is introduced similarly as in [14], i.e. by starting from the approach 
described in [10]. That is, we assume that the medium has a size L, and that the collisional interaction has to occur 
inside the medium. 
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FIG. 1: Feynman diagrams Mi,o,2,c, Afi, 0,2,1, and Mi,o,2,fl contribute to the radiative energy loss to the first order in the 
opacity. On each panel, left (right) gray ellipse represent the source J, which at time xo (xs) produces an on-shell jet with 
momentum p2 (ps). The large dashed circles ("blob") represent efltective HTL gluon propagators [15]. A cut gluon propagator 
with momentum k and color c corresponds to the radiated gluon. A cut gluon propagator with momentum q and color a 
corresponds to a collisional interaction with a parton in the medium. Specific time points are represented by Xi. The diagrams 
are calculated in light cone coordinate system, and X2 — Xq < 2L on the first and the third panel, as well as — < 2L on 
the first and the second panel, represent the condition that the distance between collisional interaction and jet production has 
to be smaller than the size L of the medium. Left, middle and right panels present three possible cuts (central, left and right, 
respectively) of the same 2-HTL Feynman diagram, all of which contribute to the 1^* order in opacity radiative energy loss. 

As in [8] , we describe the medium by a thermalized quark-gluon plasma at temperature T and zero baryon density, 
with Uf effective massless quark flavors in equilibrium with the gluons. The formalism for computing the energy 
loss in finite size dynamical QCD medium is presented in Appendices C-N, and the diagrams are evaluated in finite 
temperature field theory [16, 17], using HTL resumed propagators [17] for all gluons. To outline the calculations, 
in Fig. 1 we show three typical diagrams that have to be computed. The Feynman diagram in left panel of Fig. 1 
represents the source J, which at time xq produces an on-shell jet with momentum p2, and subsequently radiates a 
gluon with momentum k = {uj, kz, k) and exchanges a virtual gluon of momentum q = (qo, 9z, q) with a parton in the 
medium. The quark jet emerges with (measured) momentum p = {E,p2,p). We assume, as in [18], that J changes 
slowly with jet momentum, i.e. that J{p + k + q) « J{p + fc) « J{p)- To incorporate the effect of finite size QCD 
medium, we assume that the distance between the jet production and collisional interaction has to be smaller than 
the size of the medium. 

Since the produced jet can be off-shell, the Feynman diagrams shown in the center and right panels of Fig. 1 also 
contribute to the L** order in opacity radiative energy loss. They are complex conjugates to each other and they 
present the terms which will interfere with the diagram shown in the left panel of Fig. 1 , and consequently lead to the 
appearance of LPM effect (in the case of high energy jets), after all relevant contributions are taken into account. In 
addition to the three diagrams shown in Fig. 1, there are 21 more diagrams that contribute to the 1^* order in opacity 
radiative energy loss, and their calculation is presented in the Appendices C-N. Note that the calculation presented 
in this paper differs from Ref. [13] by the use of HTL gluon propagators to describe the interaction of the quark with 
the medium. The difference from Ref. [8] is that in this work we allow the jet to be on- or off-shell and the vertices 
that correspond to gluon exchange are restricted to be located inside the medium. 

Since the exchanged gluon momentum is space-like [14, 19, 20]), only the Landau damping contribution {qo < |q|) 
to the cut HTL effective gluon propagator D{q) needs to be taken into account [14, 19, 21]. The radiated gluon has 
timclikc momentum k = (w, k), so only the quasi-particle contribution at w > |k| from the cut gluon propagator D{k) 
contributes [15-17]. Since our focus is on heavy quarks with mass M ^ gT, we neglect the thermal shift of the heavy 
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quark mass. 

The effective gluon propagator lias botli transverse and longitudinal contributions [22-28]. The 1-HTL gluon 
propagator has the form 



(2.1) 



where I = (Iq, 1) is the 4- momentum of the gluon and P^ij{l) and Q^u{l) are the transverse and longitudinal projectors, 
respectively. The transverse and longitudinal HTL gluon self energies YIt and TLl are given by 
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where y = lo/\l\ and n = gTy^NJJTNjJd is the Debye screenmg mass. 

While the results obtained in this paper are gauge invariant [20], the calculation is for simplicity presented in 
Coulomb gauge. In this gauge the only nonzero terms in the transverse and longitudinal projectors are 
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As in [13, 18, 29-34], we use the same kinematic approximations, i.e. we assume validity of the soft gluon {oj <C E) 
and soft rcscattering {lu ^ |fc| ^ |q| ^ qo,qz) approximations (see Appendix A for details). In Appendices C-N we 
compute all the diagrams that contribute to the first order in opacity radiative energy loss. Once the diagrams are 
calculated, the interaction rate is given by: 



T{E) = -LMtot - -^(A/i,o + Mi,i + Mi,2), 
-'Vj iVj 



(2.4) 



here Mi^q, Mi^i and Mi_2 present the sum of all contributions in which zero, one or two (respectively) ends of the 
exchanged gluon q are attached to the radiated gluon k. Furthermore, Nj is an integrated invariant distribution of 
jets, created by the effective jet source current, and given by [18] (note that 3 accounts for the jet colors) 
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Equations (C16), (D13), (D29), (E7), (F12) and (G7) give final results for the Feynman diagrams contributing to 
Afi.Q. After adding these expressions Mi^ becomes 



D, 



{2TTf2E 



\J{P)\ 



{2ttY2E 
2 Cuas L 



\J{pt 



d^k (fq 



(27r)32cj (27r)2 q^{q^+fi^) {k^+x^ 



TT X 



dyn 



fc2 



dx d'^k d'^q 



sin(eL) 



(2.6) 



are defined by Eq. (A8) in appendix A. 



where x = M^x^ + rri^ and ^ ; 

Here and below we used [ta,tc] [tc,ta] = C2{G)CiiDii (with C2(G) = 3, Cr — |, and = 3) and defined a 
"dynamical mean free path" (see [8]) Adyn through 



dyn 



(2.7) 



where as = |^ is coupling constant, and we assumed constant coupling g. 

Equations (Hll), (18), (JIO) and (KIO) give final results for the Feynman diagrams contributing to Mi,i. After 
adding these expressions A/i.i becomes 



MiA = Di 
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X TT TT q'^{q^ + IJ? 



2k-{k+q) / sin(Ci)\ 



ik^+x)iik+qY+x) 



(2.8) 



where C = is defined by Eq. (A8) in appendix A. 
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Equations (L7) and (M7) give final results for the Feynman diagrams contributing to Afi,2- After adding these 
expressions Mi^2 becomes 
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X n TT q^{q^ + ^i^) 
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By using Eqs. (2.5)-(2.9) the interaction rate (Eq. (2.4)) reduces to 
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(2.10) 



The heavy quark radiative energy loss per unit length is obtained from the above expression for the interaction rate 
by weighting it with the lost energy uj + q{). In the soft rcscattering approximation uj -\- ~ lo, leading to: 



This finally leads to 
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(2.12) 



It is important to note that, similarly to [8] , each individual diagram that contributes to the energy loss in a finite size 
dynamical QCD medium diverges logarithmically in the limit of zero transverse momentum of the exchanged gluon, 
q^O. The reason for this divergence is that in a dynamical QCD medium both transverse and longitudinal gluon 
exchange contribute to the radiative energy loss [35] . While Debye screening makes the longitudinal gluon exchange 
infrared finite, transverse gluon exchange causes a well-known logarithmic singularity [17] due to the absence of a 
magnetic screening [36]. Remarkably, we see from Eq. (2.12) that, when the contributions from all diagrams are taken 
into account, the infrared divergences cancel, naturally regulating the energy loss rate. 

The analytical expression for the energy loss in dynamical medium that we obtained, can now be directly compared 
with the equivalent expression for static QCD medium. We below make this comparison in order to study the 
importance of the dynamical QCD medium effects. To do that, we here rewrite the expression for the DGLV first 
order in opacity radiative energy loss in static QCD medium, which is obtained in [13]: 

A£;,tat Cna, L f^^fkfq 



E TT X. 



'Stat 




with [8, 37] 



^ ^ {k+qf + x ' ^ "- ' "'"'^ ' '-2 , , ) ( 1 ik4^n\^-i-^ , I ' (2-13) 



1 1 202 1 + ^ 1 



Astat Ag Xq TT^ 1 + Ad; 



where c{nf) = 6 ^'^2 ^ i+"j/6 ^ slowly increasing function of ny that varies between c(0) ~ 0.73 and c(oo) w 1.09. 

For a typical value Uf = 2.5 (which we use in this paper), c(2.5) « 0.84. 

We see that, similarly to the case of infinite QCD medium [8], Eqs. (2.12) and (2.13) arc remarkably similar, up to 
two important differences: The first is an 0(15%) decrease in the effective mean free path 



(2.15) 
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which increases the energy loss rate in the dynamical medium by 0(20%). The second difference is a change in the 
effective crossection, which in the energy loss rate is reflected by the replacement 
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(2.16) 



As discussed in [9], these differences lead to a significant increase of the heavy quark energy loss rate in dynamical 
compared to static QCD medium. 

By using the above results, we can now discuss the following two issues: 1) Comparison between dynamical and 
static energy loss results, 2) Comparison between energy loss results in Bethe-Heitler limit and finite size QCD 
medium. Regarding the first comparison, one should note a remarkably simple mapping between dynamical and 
static QCD medium. That is, the expression for energy loss in dynamical QCD medium can be obtained from the 
expression for the energy loss in static QCD medium by simply replacing the effective mean free path and the effective 
crossection by the appropriate expressions given above. The simplicity of this substitution rule is surprising, given 
the complexity of the calculations and their different structure for static [13] and dynamical media. In particular, 
one should note the infrared divergences in the dynamical case which cancel only after summing all the diagrams 
from Appendices C-N, but don't arise at all in the static case. Regarding the second comparison, the expressions for 
both dynamical and static energy loss in finite size QCD medium are significantly different from the corresponding 
expressions in Bethe-Heitler limit [8]. However, despite this difference, the same simple substitution rule is found to 
apply, suggesting a possibly general mapping between static and dynamic QCD media. 

We also note that the study presented here considers a finite, optically thin dynamical QCD medium (QGP), 
extending the DGLV approach [13, 18] to include parton recoil. In this sense it is complementary to the work by 
Arnold, Moore and Yaffe [34] who study energy loss in an infinite, optically thick QGP. We note that the AMY 
approach [34] yields the same form (2.16) for the effective cross section in a dynamical QCD medium as found here 
(see also [38]), further supporting our conjecture above. 



III. INCOHERENT LIMIT 



First of the two relevant limits of the Eq. (2.12), which we consider in this paper is the incoherent, short formation 
time limit. This limit is relevant for heavy quarks at lower jet energy regions, and it can be extracted from the 
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where the second step is obtained after angular integration. Furthermore, under the assumption that as is not running, 
Eq. (3.1) can be further analytically integrated over < ]fc] < fcmax where k^ax = 2E^/x{T^~x) [13]. We obtain 



AE, 



dyn 



E 



CrUs L 

7r2 Ady: 



— ^ [ dxd'^qJdyn{q,x) , 

Adyn J 



(3.2) 



where 



q2 (qt2 _|_ ^2-) 
X 



q 



v/q4 + 2q^x-kl,.) + (^Lx+x) 



In 



A'2 

max 



g2 + ^qi + 2q-^{x-kl,^) + {kl,^+x) 



2 (X + A;^ax) 



{q" + 2x) 



In 



c+X 



(qn3x) + ,/l+|$(g^+x) 



_4x 



(q'-fc^ax+3x) + V?" + 292(x-fc,LJ + {kla^+xY 



6 



By comparing Eqs. (3.1)-(3.3) with the equivalent expressions from [8] (see Eqs. (2.9)-(2.12)), we see that, though 
similar, the corresponding expressions are not the same. The reason is that, in [8], we considered the Bethe-Heitler 
limit, which considers on-shell quark jets produced at remote past, while in our study, wc allow that quark jets can 
be both on- and off-shell. 



IV. LPM LIMIT 



The second relevant limit that we consider is the Non-Abelian analog of the Landau-Pomeanchuk-Migdal limit. 
This limit is relevant for highly energetic jets, where destructive interference effects reduce the energy loss relative 
to the incoherent limit. The limit is obtained from the Eq. (2.12) when « 2E oo. In such a limit finite mass 
effects arc negligible. Additionally, /cmax ~ 2EyG^(T^-x) — > cx) as well, which enables us to introduce a substitution 
fe' = fc + q in Eq. (2.12). With these simplifications, Eq. (2.12) becomes 

A^dyn ^ Cn^_L_ f d-'k' d?q fi^ 2q-{q-k') ( 

E IT Adyni ^ TT q^[q^+ti^) k''-{k' -qf \ MLL j 

^ a. r r'^^ V .(1.1 H^'i) fi - , (4.1) 



71- Adynio Jo Jo q^{q^ + i-L^) \ ^ 



2xE 

where in the second step we performed angular integration. To proceed further, we observe that the derivative over 

dL' 



distance L of the above expression (i.e. 4^) is equal to 



dyn _ 2-^^- / da; / / dq^^—^ ^ 6l(|q| - |fe'|) 1 - cos(— — ) 



E dL TT Adyn 7o Jo Jo q^(g^+A*^) V 2x£' 

' (. - c. m . m)-] . (4.2) 



TT AdynJo Jo 9^(9^ V \2Ex J \2Ex 

where in the second step we performed integral over and 7 « 0.577216 is Euler's constant and Ci{y) gives the 
cosine integral function. 

Finally after performing the integration over x, we obtain 



E dL TT AdynJo q^{q^^lJ?) \ \2E ) 8E 

where 

Ziy) ^ 7 - J cos(y) - C^{y) + ln(y) + - 0. (4.4) 

2 Zy 2 y^Q 

Therefore, for the asymptotically large jet energies, the Eq. (4.3) reduces to 

E dL 2E Adynio g2+^,2 2E Ady„ ' V A^" / 

where we used gmax = ViET [39] . Finally, ^ then becomes 

E AE Adyn / ■ 

From the above expressions, we see that at asymptotically large jet energies we obtain quadratic thickness dependence 
for the energy loss, that is we recover LPM limit. Therefore, for highly energetic jet, finite size corrections implemented 
in the calculation presented in this paper simulate the destructive effects of LPM interference in an infinite QCD 
medium [5]. This behavior is expected [18] since the nuclear medium has finite dimensions that may be small 
compared to the jet radiation coherence length, especially in the case of light partons or high jet energies. Due to this, 
in finite size media the basic formation time physics developed by LPM [12] leads to destructive interference effects 
on the quark quenching. 
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Furthermore, by applying the same procedure, it is straightforward to obtain that, for asymptotically large jet 
energies, the radiative energy loss in static QCD medium becomes 

so that ratio between energy losses in dynamical and static QCD medium approaches 

stat in ^ A,. 



hm -rrr^ = 1™ ir^ = t ■ (4.8) 

E^oc Ai^stat -E^oo Adyn 1 Adyn 



hcrcfore, we conclude that, at asymptotically large jet energies, approximation of the medium by a random distribution 
of static scatterers becomes valid, up to a multiplicative constant which can be renormalized. This is consistent 
with what would be expected from established BDMPS results [5]. 



V. CONCLUSION 



In this paper, we developed a theoretical formalism for the calculation of the first order in opacity radiative energy 
loss of a fast quark traveling through a finite dynamical QCD medium. We obtain a closed analytical expression 
for the energy loss in dynamical medium. The obtained result is convergent, despite the fact that each individual 
contribution to the energy loss is infrared divergent. Furthermore, the energy loss has a nontrivial dependence on the 
size of the medium, which depends on both mass and energy of the quark jet. The finite size effects are found to be 
most important in the ultra-relativistic limit and they effectively reproduce the effects of destructive LPM interference. 
Another interesting limit is an incoherent (GB) limit, which is reproduced for heavy quarks with moderately small 
jet energies. 

The study presented here considered the radiative energy loss in finite size dynamical QCD medium up to the 
first order in opacity. On the other hand, in static QCD medium, radiative energy loss up to all orders in opacity is 
obtained [13, 18]. Simphcity of the mapping between the l*** order in opacity static and dynamical energy loss, implies 
that the same mapping might be generalized to higher order in opacity as well. That is, we make a conjecture that, 
to obtain the energy loss expressions in dynamical QCD medium from the existing static QCD medium expressions, 
one (only) needs to replace effective mean free path and effective crossection from static QCD medium, with the 
corresponding expressions from dynamical QCD medium. However, to prove this conjecture is very non-trivial, which 
will be a subject of further research. 

The measurement of the heavy flavor suppression is in the current focus of intensive experimental efforts, and these 
measurements are expected to became available soon at the upcoming high- luminosity RHIC and LHC experiments. 
As already mentioned, particle suppression is a consequence of the energy loss. Our study, which incorporates 
dynamical effects in realistic finite size QCD medium, enables us to provide the most reliable computations of the 
energy loss in QGP. Our future goal is to use these energy loss calculations to make accurate theoretical predictions 
for the heavy flavor suppression. These predictions can then be directly compared with the upcoming experimental 
data, in order to test our understanding of QGP, and to further study the properties of this novel form of matter. 
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APPENDIX A: NOTATION AND ASSUMPTIONS 



In the following appendices, the calculation will mostly (with the exception of G |_ propagators), be done using the 

light cone coordinate system [40] . This coordinate system is appropriate for systems moving with almost the speed of 
light. It is obtained by choosing new spacetime coordinates [x~^ ,x~ ,x], related to the coordinates in the laboratory 
frame {t, z, x) by {x is the transverse coordinate) 



x^ = (t + z), X = (i — z). 



(Al) 
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In the same way the hght cone momentum [p'^ ,p ,p] is related to the momentmii in the laboratory frame [E,pz,p\ 
by {p is the transverse momentum) 

p+ = {E+p,), p- = {E-p,). (A2) 

In this paper we consider a quark of a finite mass M, produced inside the finite size QGP, at some initial point xq. 
The jet is assumed to have a large spatial momentum p' 3> M, and it can be both on-shell and ofF-shcU. This is in 
contrast to the Bethe-Heitler limit, which we considered in [8], where an on-shell jet was created at — oo. Further, we 
choose coordinates such that the momentum of the initial quark is along the z axis: 

p' - [E'+,p'-, 0] , (A3) 

We are interested in the radiative energy loss to first order in the opacity, so we study the case in which the quark 
exchanges (in arbitrary sequence) one virtual gluon with space-like momentum 

q = [q'^,q^,q] = (90, q) = (go, qz, q), qo < |q| (A4) 

with a parton in the medium and radiates one (medium-modified) real gluon with time-like momentum 

k^[k+,k-,k]^{ka,k) = {oj,k„k), ka>\k\ (A5) 
into the medium. The on-shell quark jet emerges with 4- momentum p^. 

Similarly to all energy loss formalisms developed so far (see e.g. [5, 13, 18, 29-34]), the calculations presented in 
Appendices C-N are developed under the assumption of a perturbative high-temperature QGP. We note that, strictly 
speaking, this assumption may not be directly applicable to the case of a strongly coupled QGP or in cold nuclear 
matter. The validity of this assumption for the QGP created in ultra-relativistic heavy-ion collisions can be tested by 
comparison with upcoming experimental data from the RHIC and LHC heavy-ion programs. 

For the computation of the Fcynman diagrams given in Appendices B-D we will need G^^(.t), G'I_(.t) and G |-(.t) 

propagators for the quark jet p, the radiated gluon k, and the exchanged gluon q. These functions are derived in 
Appendix B. 

As in [13, 18, 29-33], we assume validity of the soft gluon {uj<^E) and soft rescattering (jql ~ |A;| <C fc3) approx- 
imations. These approximations become more reliable as the temperature increases and are expected to hold well 
at the LHC, while at RHIC their application should be further validated. Together with conservation of energy and 
momentum {p' = p + k + q) they yield 



k ^ 



E^,P = , p 



(A6) 



P : 

In [8] we showed that it is reasonable to assume that q^ has the same order of magnitude as |<7|. Since |fc| ^ k^ and 

E+ ■ 



qz ~ \q\ ^ \k\, we then also have qz <^ kz- Thus k^+qz ~ kz and Pz+kz+qz ~ Pz+kz ~ Pz+qz ~ Pz- Defining x = 



we can therefore also assume that 

_ k+ {k + qy 



E+ {E + q)+ ' 
By using x, we further define C ^i^d Q in the following way 



(A7) 
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xE+ 
{k+qf 



xE^ 



Finally, by using that ^,C^ Ifcl'^^l?!, ^-^d that ~ |<7| ~ |fe|, we obtain X ^ qz leading to 

q- +£. = q^-iqz-O-q^'-qz 

q-+C = q° -{qz-C)^q"~ qz 

g-±(C-0 = q°-iqz±iC-0)^q"-qz- (A9) 
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APPENDIX B: DERIVATION OF THE PROPAGATORS IN LIGHT CONE COORDINATE SYSTEM 

In this appendix we present in some detail the derivation of the propagators G±± in light cone coordinate system, 
which are needed for the calculations presented in the following Appendices. To derive these functions we start from 
Kallen-Lehmann theorem [41]: 

^ ^ i°° - / - ■ 

where a, b can be + or — , and Pabis, ^^'^ Pab{si ^^'^ spectral functions. We can now decompose the propagator as 

Gabix) = eix+)Gt,{x) + ei-x+)G-,ix) (B2) 

Note that 

l-x = ^l^x^ + ^l^x+ - Ix (B3) 
f = l+l- ~f , (B4) 

We first concentrate on G^^{x). To obtain G^f^{x) we take poles from Eq. (Bl). Pole P — s + ie = contribute to 
G^j(x-) only if /+ > 0, while pole — s — je = contribute to G^j(x) only if < 0, leading to 

Gtbi^) = j J2ny^'^'' j, dsS{P~s){9{l+)p,,{s,l)-e{-l+)p,,{s,l))^ j j^e-^^-^Gt,{l), (B5) 

where 

G+b(0 = e{l+)pab{l\l) - e{~l + ) pabif, I) , (B6) 

where 0{l^) and 6{~l^) are unit step functions. 

To proceed further, we use the following expressions 

p++(/^^) = -p^^{f,i) = {e{f) + .f{f))p{i) 
P-{i\i) = -p++{iM)^{e{-f) + l{i'))p{i) 
P+-{i\i) = -p+^{i\i) = .f{f)p{i) 

p-+{l\l) = -p^+{l\l) = {\ + f{f))p{l), (B7) 

which can be straightforwardly derived from [41]. Here 9{f') is unit step function, /(/") = (e'"/-'"— 1)^^, where T is 
the temperature of the medium, and p{l) is spectral function. 
Then, by using Eqs. (B6) and (B7), we obtain 

G++(o = -e{i+) p++{i) + e{-i+) p__(0 = {6{i+) 6{f ) + e{-i+) e{-f) + /(?")) pH) 

Gt^il) = -6(1+) p__(0 + 9i-l+) p++{l) = {9(1+) 9i-f ) + 9i-l+) 9if) + /(?«)) pH) 
G+_il) = -9{l+) p+41) +0{-l+) p+-{l) = P+-{1) = fil"") P{1) 

Gt+{1) = -9{l+) p^+{l) + 9{-l+) p^+{l) = p_+(0 = (1 + /(?")) p{l) . (B8) 
Furthermore, it is straightforward to show that 

G++(0 = g:_(0 

G:+(0 = G+_{1) (B9) 

In this paper we will need only G\^{x), Gl„(x) and G ^-(a;) propagators for the quark jet p, the radiated gluon 

fc, and the exchanged gluon q. We will first derive these functions for the exchanged gluon q. To do this, lets first 
note that, as noted in the previous appendix, in this paper we assume high temperature QGP and soft-gluon, soft 
rescattering approximation. In such a limit, for the exchanged gluon g, /(q°) « ^ ^ 1, reducing the Eqs. (B8) 
and (B9) to 

G++{q) « Gl_iq) « G+_{q) « Gl+{q) « G:_(<z) « GZ+{q) « /(q°) p{q) . (BIO) 
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By using Eq. (BIO), it is then straightforward to obtain that corresponding propagators for the exchanged gluon are 
given by 

(x, - X,) = DZ't'ix, - X,) « Dl^_{x, - X,) « j ^D>iq)e-"^(^'-^^^ , (Bll) 
where D^[q) is the effective 1-HTL cut gluon propagator for the exchanged ghion [8] 

0>M) . .(1 - |) (1 + /(,„)) 2I.„ (-,2^ + Ji^l ^ (B12, 



Here ^t{q) and ^hiq) (see Eq. (2.2)) are transverse and longitudinal self energies. Note that exchanged gluons are 
spacelike (see Eq.(A4)). Therefore, in the cut 1-HTL exchanged gluon propagator, only Landau damping contribution 
from the gluon spectral function contribute to the above expression (for more details, see [8]). 

We will now concentrate on the propagators for the radiated gluon and quarks jet. For radiated gluon and quark 
jet /(Zo) = (e'o/^-l)-! < 1, reducing the Eq. (B8) to 

Gi+(0 = {e{i+)9{f) + e{-i+)e{-f))p{i) 

G+_il) - 

Gt+(/) - p{l) (B13) 

To proceed further with propagators for the radiated gluon propagator, we note that, for the radiated gluon with 
momentum k, the longitudinal contribution can be neglected relative to the transverse one. Also, for the transverse 
gluon the self energy Urik) can be approximated by m^, where nig « /x/v^ is the asymptotic mass (sec [15]). These 
approximations are true in the soft rescattering limit u 3> |q|^|fe|~(7T which we use in this paper. With these 
approximations the HTL gluon propagator for the emitted gluon can be simplified to [15] 

D,.{k)^-^TA^^^, (B14) 

where P^i^ is the transverse projector. 

Then, by using Eqs. (B13) and (B14), relevant radiative gluon propagators become (see also [40]): 

/rlh+rPk 
^^e(k+)P^%k)e-^^^^^-^^) (B15) 

___pA<.(fc),-./c(.,-.,) (B16) 



where w « y + m?g and fc^ = + /c^ . 

Note that we here assume that gluon mass is given by the expression mg w /i/-\/2 regardless whether the gluon 
is radiated inside or outside of the medium. Strictly speaking, including the finite size effects on the radiated gluon 
would assume that the above gluon mass is valid only for gluons radiated inside the medium, while for the gluons 
radiated outside the medium the gluon mass should be equal to zero (see [33]). However, based on [33], neglecting 
the finite size effects on the radiated gluon is expected to be a reasonable approximation, and we adopt it here to 
simplify the calculations. 

Finally, we now concentrate on the relevant propagators for the quark jet. By repeating similar procedure as for 
the radiated gluon, we obtain 

Ai+(^. - X,) - A:_(x. ~ X,) ^ j ^g^0(p+)e-P(-.-.) (B17) 
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APPENDIX C: COMPUTATION OF DIAGRAM Mi^o^.c 

In appendices (C) - (G) we present in some detail the calculation of the diagrams where both ends of the exchanged 
gluon q are attached to the heavy quark, i.e. none is attached to the radiated gluon k and no 3-gluon vertex is 
involved. In this appendix, we start with the calculation of the diagram shown in Fig. 2. 

Here and later the diagrams are labeled as follows: In Mi^ij^c, 1 denotes that these diagrams contribute to the 
energy loss to first order in opacity; i denotes how many ends of the virtual gluon q are attached to the radiated gluon 
k; and j labels the specific diagram in that class. Letter C denotes that we consider central cut of the diagram. In 
the next chapters, letter R [L) will denote that we consider right (left) cut of the Feynman diagram. 




FIG. 2: Feynman diagram Afi,o,i,c, contributing to the radiative energy loss to first order in opacity, labeled in the same way 
as Fig. 1. 

1. We will first calculate the cut diagram Afi.o.i.c: 

X A„+(a;3 - X2)v;^{x3) AI_(.T4 - X3)v:^{x4) Al_(x5 - a;4) J(a;5) 

X e{x+ - x+) 9{x+ - x+) e{2L - {x+ - x+)) e{x+ - x+)e{x+ - x+) e{2L - {x+ - x+)) , (CI) 

where J is the source of a jet, A correspond to the jet propagator, D to gluon propagators and v to vertices. By 
using expressions for the propagators from the previous section, we obtain: 



/ '^ poo f 

\{d^x,J{x^) / / 



dp2d^P2 
(2^)3 2p+' 



-ip2-ixi-xo) 



-iq{x4 —xi) 



dpjd^pi 



X (-1) 



(2^)3 2p^ 
d^k 



(27r)32w 



Ppa{k) 



^ — ik-{x3~X2) 



d^p 



°° dp+d^p. 



oo poo 



(2^)3 2p+ 

^2 ^ 

dpfd^Pi f d'^k 



X J{x5) e{x+ - x+)9{x+ - x+) e{2L - - x+)) e{x+ - x+) 6i(x+ - x+) e{2L ~ (x+ - x+)) 



(2^)3 2p+ 

+ 

3 



n 



d^p 



dS 4 



fJii2TT)3 2p+ J {2ttY2uj J {2TTf2E J {2tt) 

X [pi+P2rDiM){p^+Pir {p+piYPp.{k){p+p^r I 



g tatctc^a 



(C2) 
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where 
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I ^ Yi d^^^^i^t - ^o) ^("^t - ^t) f (2i - (2;^ - x^)) Oixi - x+) 9{x+ - xj) 9{2L - {xj - x+)) 

i=0 

X ^-■i.P2-{xi-xo) ^-iq-(xi-xi) ^-ipi-(x2-xi) ^-i{p+k)-(x3-X2) ^-ips-ix^-x^) ^-ipi-{x^^-Xi) J^j;^) J{x^) 

= I J(p)|2 (2^)3 5{{p2 -p-k- q) + )5\p2 -p^k q) {2nf S{{p, - p - k) + )S^{p, - p - k) 

X (27r)3 ^((p3 -p- k)+)S\p3 p k) {2nf 5{{pi-p-k- q)+)5^{pi ~ p ~ k - q) h , (C3) 



where 



Jo Jo Jo 



Here we defined x[ = xi — Xo, x'2 = X2 — Xi, x'^ = X3 — Xi, x'^ = X4 — X5. 
By applying S functions from Eq. (C3), and by using 



Pi 



pj + AP 



k ^ (C5) 



we obtain (note p + k+q = ^ p + fc = — q) 



Then m the soft gloun, soft rescattermg hmit we obtain (note x = and that ^ and are defined in Eq. (A8)): 

{pi-p- k)~ = {p3 -p-k)~ = = 

{P2-p-k) = ^^^^^ = -C (C7) 



leading to 



^6 ^ sin'((p2-p-fc-g)" f) 

{{p.2-p-k-q)-f 



h = „^ ,.^,. 4 -;^"- " - (C8) 



In the paper [14], it was shown that finite size effects on collisions are negligible, i.e. for coUisional parts 

^ sm^((p2-p-fc-g) f) ^ -p-k-q)-)^2-KL S{q- + C) « 27tL <5(q" - g,) (C9) 

[[p2-p-k-q) ) 

where in the last step we used soft gloun, soft rescattering approximation, i.e. Eq. (A9). 
By using Eq. (C9), Eq. (C8) reduces to 

I^^Jj^^2^-L5{q'-^^)- (CIO) 
Similarly as in [8], for highly energetic jets 

{p+PiYPpAk){p+P3r « -4^ 

{pi+P2rP^.M){pz+Pir ~ -{Pi+P2YQ^u{q){P3+Pir ~ ^ (cn) 
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By using Eq. (Cll) and Eqs. (B12), {pi+P2)''D>^{q){p3+p4Y becomes 



{pi+P2rD>MiP3+P4r ^ bi+P2r2im(-^%^ + #^)(P3+P4r 



(l-^)/(go)i?+%2Im 



,2 / 1 



1 



W-IiLiq) q'-nTiq)J ^^^^^ 

Finally, by using Eqs. (C3), (C4), (CIO), (Cll) and (C12), and after performing integrations over pi, p2, P3 and 
P4, Eq. (C2) reduces to 



d^p I J, d^k 

(27r)32£;' (27r)32w (fc2 + ^) 



4, 



2 ^9 



where 



—--j2TT6{qo - g^) /(go) ^ 2Im 



q'-ULiq) q'-IlT{q) 



In high temperature limit Ig reduces to (see Appendix C in [8]): 



T 



d^q fi^ 



(27r)2 q^q^+fi^'' 



Finally, Eq. (C13) becomes 



'3- '• d^k d^q 



d^p 



{2TTf2E 



\J{p)Y 



(271)320- (27r)2 q2(q2+^2) (fc2+^): 



(C13) 
(C14) 

(C15) 

(C16) 



APPENDIX D: COMPUTATION OF DIAGRAMS Mi,o,2,c, Mi,o,2,h AND Mi,o,2,i 
In this appendix we present in some detail the calculation of the diagrams shown in Fig. 3. 



1^1,0,2,0 




J(X5) J(Xo) 



(Hb: 



J(>f5) 




FIG. 3: Feynman diagrams Mi,o,2,c, M\fi^2,L and Mi,o,2,fl contributing to tlie radiative energy loss to first order in opacity, 
labeled in the same way as Fig. 1. The same figure is also presented as Fig. 1 in the main text, and is repeated here for 
completeness. 

1. We will first calculate the cut diagram Afi.o,2,c, shown in the left panel of Fig. 3: 



Mi.o. 



2,C 



/ d'^Xi J{xo) A+_|_(a;i - xq) u+(a;i) D^!'^(a;4 - xi) A|_^(a;2 - xi)v^{x2) D''\{xs, - X2) 

i=0 

X A_+(a;3 - a:2)tv (xs) Al_(x4 - X3) w^(x4) Al_(x5 - 2:4) Jix^) 

X 0{xt - x^) 9{x^ - xt) 0{2L - {xt - x^)) 9{xt - x^) e{xt ~ xj) e{2L - [x^ - x^)) 



n 



dpfd'^Pi f d^k f d^p f d'^q 4 



-00 JO (27r)3 2p+ 7 (2^)32^7 (2^)32i? 7 (2^)4 
X {pi+P2rP^.{k){p3+P4r [p + piYD>M{P+P5Y I 



g tototot^ 



(Dl) 
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where 

/ 



5 

/ n d'^x^eixt - 4) e{x+ - x+) e{2L - (x+ - 4)) - 4) 0(x3+ - x+) 0(2l - (4 - 4)) 

X Q-iP2-ixi-xo) ^-ik-{xi-xi) ^-~ipi-{x2~xi) g-i(p+g) ■ (x3 -3:2) g-ip3-(a:4-2;3) ^-ipi-ix^-Xi) J(^r^^'^ J{x^) 

I J(p)P (27r)3 5{{p^ -p,- k)+)d^p2 -Pi-k) {2nf d{{p^ -p~ q)+)5^{pi p - q) 



X 



(27r)^ 5((P3 - P - g)+)<5'(P3 - P - g) (2^)^ <5((p4 - P3 - fc)+)<^'(P4 -p3-k)h, (D2) 



and where 



2L rx'+ r2L 



Jo 

Here we defined x[ ~ xi — xq, X2 ~ X2 — xq, x'^ = x^ ^2:5, x'^ ^ X4 — x^. 

By applying 6 functions from Eq. (D2). and by using Eq. (C5), we obtain Pi — p^ and P2 = 
Ii then becomes 

15 g-*(P2-p-fc-g)"i _ 2^ / pi{P2-p-k-qy L _ ^ pi{pi-p~q)^ L _ ^\ 

{P2 - Pi - k)-)^ ^ {p2-p-k-q)- \{p2-p-k-q}- [pi - p - q)- j 

1} 



(pi-p-g) y{p2-p-k-q) [pi-p-qY 
16 ,. sin^((p2-p-fc-g)'-|) A {p2-p-k-q)- e'^P^'P-iy ^ - 1 



{4 



sin^ 



((pi - / _ e'ipi-p-i)-L _i {p2-p-k-q)- \ 

\pi-p-q)-f V iPi-P-iy e'iP^-P-k-D-L 



{{Pi -p-> 

In the paper [14], it was shown that finite size effect are smah for coUisional energy loss, i.e. for coUisional parts 
the following approximation can be used 



sm {{p2-p-k-q) 2) _ „ .c.. , 

4 — -7^ ; , ,2 ~ 2TTLd{{p2 - p - k - q) ) 

[\P2-p-k-q) ) 

^^^-^-^-^^ ^2.LSiip,-p-q)-) (D5) 

((pi-p-q) ) 

By using Eq. (D5), Eq. (D4) reduces to 

327TL j , , e^(pi+fe-P2)^i _ l\ / e-^(pi+fc-P2)-L _ 1' 

327rL 



(pi + fc - P2) )^ 



P2 



p-k-q) ) + 5{{p,-p-q) j 

(pi + fc -p2)^L) 
{pi + k - P2)^L 



In soft gluon, soft rescattering approximation 

iPi^P) « 

(Pi+fc-p2)- (D7) 

Note that {pi — p)^ ^ ^ <C |fc|, |q|,g2, leading to 

{5[[p2-p-k-q)-)+5{{p^-p-q)-)) « 25{{p2- p-k-q)- )^25[q^-q,), 
{K{P2-p-k-q)-)-5{{pi-p-q)-)) ^ 6{qa-q,+i)-5{qo-q,). (D8) 
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(D9) 



By using Eq. (D8), Eq.(D6) finally reduces to 

^ ^ 16x^E+^ r^^, , A sm{£.L)\ .( „ , \ cos(CL) - 11 
^1 « (fe2 + ^)2 |2<5(<?" - g.) (^1 - j +4^(90 - 9. + - '5(90 - 9^)) 1 

Similarly as in previous section, for highly energetic jets and in high temperature limit 
By using Eqs. (D2), (D9) and (DIO), and after performing integrations over pi, p2, Vz and p4, Eq. (Dl) becomes 



(27r)32w (fc2+x)2 
sin(eL)\ . cos(gL) -l 



X ^2(1-^^)/, + . -^^-/ - J, , (Dll) 



where Iq is given by Eq. (C15), and 



In the last step we used that ^^^-^^j^^^^^ |qo=g, is an odd function of Qz- 
Finally, by using Eqs. (C15) and"(D12), Eq. (Dll) reduces to 

. f d^p , , ,,9 f (fk cPq ( sin(£L)\ , 

2. We will now calculate the cut diagrams Mi^o,2.l and Mi_o,2,R, shown in the central and right panel of Fig. (3), 
respectively. We start with Mi o,2,_r: 

X A+_^(a;3 - .T2) v+(.T3) A_+(x4 - 2:3) ^^(.14) AI_(.T5 - X4)J{xz) 
X 0(x+ - x+) 9{x+ - x+) e{x+ - x+) e{2L - {x+ ~ x+)) e{x+ ~ x+) 



OC pOC 



n dptffpt f <Pk f d^p f d^q 4 , , , , 



fJi (2^)3 2p+ 7 (2^)32a; 7 (27r)32i; J (2^)^ 
X (P2 +P3)^/^M-(fc)(p+P4)'^ {Pi +P2YD>{q){P + PiY I . (DM) 



where 



5 

/ = / n d'a;,0(4 - 4) 0(4 - 4) e{x+ - 4) 0(2i - {x+ - x+)) 6{x+ - 4) 

1=0 

X g-ip3-(2:i-a;o) ^—ik-(xi-xi) ^—ip2-(x2-xi) ^-i{pi+q)-(x3-X2) ^—ip-{xi-X3) ^-ipi-(x^^-Xi) J (^^^^ Jix-^) 

= \J{p)\^ {27TfS{{p3-p2 - k)+)5^p3-p2 - k) {2nf5{{p2-p)+)5\p2-p) 

X (2^)3 5{{pi ~p + q)+)S^{pi -p + q) (27r)3 ^((^4 - p - k)+)6\p4 - p - k) h , (D15) 

and where 

p2L c^''! 

h= d4+e-*(P^-P)"^^+ / dx'+e-i'^P^-P^-'''^~<^ / dx'+e-i'^P^-P+'^^''''^^ / dx^+e5(P''-P-'=)"<'' (D16) 
Jo Jo Jo Jo 

Here we defined x'^ — xi — xq, x'2 = X2 — xq, x'^ = x^ — X2, x'^ = xa ~ x^. 



By applying 5 functions from Eq. (D15), and by using Eq. (C5), we obtain 

P2 ^ ip2 -py = 0, 

{pi -p + qy q^. 

Note that in the second equation we used p = k and (p + fc)+ w . 
Ii then becomes 



JQ Jo Jo Jo 



By using Eqs. (D18), (D15) and (D14) il/i,o,2,_R becomes {y = ^) 

(Pp , , .,2 f d^f^ 



Ml^0,2,R = -ALgUctatatc / ,^ ^ |J(p)|^ 



sin^L .1 — cos^_L\ I' d^q 

^ ) J W) 



X (i_^„,±^±^) I ^yp + p,rD>Aq)ip + P.r I dye-^^-y 



where, as in the previous sections, we used (p+Ps)'^ Pi_iu{k){p+p3)'' 
In the same way, it can be obtained that Mi q 2,l is equal to 



)32a> (k^ + x) 



2 



( sin£L 1 — cosfL\ ( d'^q , , w f°° . 



Mifi^2,R + A^i,o,2,L then becomes 

d^p , J, ^,2 f 



d'^q 



(27r)32w {k^+xY 
j{p + PiYD>,{q){P + Pirh, 



where 



(2^) 



sin^L l-cos£L\ [ sin^L l-cosfL\ , 

sin^L\ „ , „ l-cos^i f°° , . ^ ^ 
1 - — ^ 27r5(g-) - 2 — ^ / dy sin q^y . 



Similarly as in previous section, for highly energetic jets 



(p+p,)'oi(<rt(p+p,)" » «(i-|)/(*)i=«|;2i,„(-^--^) 

Finally, by using Eqs.(D22) and (D23), Eq. (D21) becomes 

A f d^p ,,o /■ S'k fe^ / sinfL 



(27r)32i; ' '"^^ 7 (27r)32w (A;2 + x)H 
4 1 d3p , ^, ,,2 f d^k fc2 1 — cos^L 
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where Iq is given by Eq. (C15), and Jq is given by 



J„ 



dz 



dqo dq^ d'^q 



{sin{qoz) cos(g^z) - cos{qoz) sin((j^z)) F{qo,q^, q), 



JO J (2^) 
In the above equation, we defined F{qQ, q^, q) as 

r / 2Imn^(g) 2IninT(g) 

Wo,gz,gj- .2^^ U'z'-Rcni(g))2 + imnL(q)2 (g2 _RenT(g))2 + imnT(<z)^ 

^'(g'Oi ^) is even function of both go and g^, leading to 

/ d(7osin(go^)^'(go,9z,g) = 



/CX3 
dq^sm{qzz)F{qo,q^,q) = 0, 
-oo 



which consequently leads to 



J, = 0. 



Finally, by using Eqs. (C15) and (D28), Eq. (D24) becomes 

(2^ "•^(^^l' J (2^ (2^ g2(g2+^2) 
fc2 / sin(CL) \ 



(D25) 



(D26) 



(D27) 



(D28) 



(D29) 



APPENDIX E: COMPUTATION OF DIAGRAMS Mi,o,3,c AND Mi,o,4,c 

We will now calculate cut diagrams Mi_o,3,c and Afi,o,4,Cj shown in the Fig. (4). We start with Mi_o,3,c: 

I I 



,0,3,C 



"^1,0,4,0 




;f|-xJ<2L 




xJ-x^<2L 



FIG. 4: Feynman diagrams A'fi,o,3,c and Mi,o,4,c contributing to the radiative energy loss to first order in opacity, labeled in 
the same way as Fig. 1. 



i=0 

X A_+(a;3 - X2)v^{x3) Al_(x4 - xs)!^ (.T4) Al„(a;5 - 2:4) J(x5) 

X e{xt - x^) 0{xt - x+) e{2L - (.Ti - xa)+) e{xt - xf) 0{xt - xt) 0{2L - (xg - 0:5) + ) 

_ _ r ruMfp^ ( f f ^'g 4, , , , 

looio £|(27r)3 2p+y (2^)32c.y (2^)32i? 7 (2^)^ ^ 

X bl +P2)"i?>,(g)b + P3)'' (p + piYPp^ik) (p3 +P4)"/ , (El) 
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where 



I = [Y[ d*x,e{xi - x^)0{xt - xi) e{2L - {xi - xo)+) e{xt - xj) 9{x^ - xj) 9{2L - {xs - xs)^) 

X g-iP2-{'J:i-xo) ^-iq-{xa-xi) ^-ipi-{x2~xi) ^-ik-(xi-X2) ^-ip-[x3-X2) ^-ip3-(xi-X3) ^-ipi-(x;,-Xi) J(^xo) J{x^) 

= I J(p)|2 (27r)3 6{{p2 ^p^k~ q)+)5Hp2 -p-k-q) {27rf 5{{pi - p - k)+)5\pi - p - k) 

X (2^)3 5((p3 -p- q)+)5\p3 p-q) (27r)3 5{{p^ ^ p., - k)+)6\pi - p^ - k) h , (E2) 

and where 

poo p2L p2L pxs 

h ^ c;a;H-g-f(pi-p-fc)-4+ / rfa.'+e-3(P2-p-fc-9)-^'i+ / d^g+e^^P^-P-'^'^a"^ / dx'+e^'^P^-P^-'''^''''^* 
Jo Jo Jo Jo 

_15 p-i{P2~p-k-q)L _ 1 ( pi(pi-p-k-q)- L _ i pi(P3 -p-<?) " i _ 1 

J J I (E3) 

{Pi - P - k)- [pi - P3 - k)- {p2- p-k-q)- \iPi- p-k-q)- {ps - P - q)^ J 

Similarly as in previous sections, by using S functions from Eqs. (E2) and (C5), we obtain p^ = P2 and 

(pi-p-k)- ^ {P4 - p3 - k)- ^ . (E4) 

Ii then becomes 

^ 16 4sin2((p2-p-fc-g)-f) f {p2 - p - k - q)- e*(P3-p-7)-L _ 1 1 

^ ^ {{P2 -p-k-q)-)^ 1 g-i(p2-p-fe-g)--L _ 1 (p^-p-q)- J 

16 , x_x r , (p-,-p-k-q)- e'{(P2-p-fc-9)-+ai _ 



--^27tL6{{p2 - p-k-q)- )[l 



^2 ^'-^ ^ U (p3_p_fe_,)-^o e-»(P2-P-*--«)"-^ - 1 {p2-p-k-q)- +^ ) 

16„ N f-, sin(^i) .cos(^L)-n ^^^^ 
■2T:LS{qo-q,) ll 7^ + 1 >■ (E5) 



where similarly to previous sections, in the second step we used Eq. (A9) and the fact that finite size effects are 
negligible for coUisional interactions. That is, we used '^^^^(^pf-p-k^qyy * 2'!rL5{q- — ^) w 2'!TL5{qo — q^)- 
Finally, by using Eqs. (E2) and (E5), Eq. (El) becomes 

M1.0.3.C = -ALTgHaUJ, [ -fjL^\J[p)\^ I' 



(27r)32£; ' J (27r)32w (k^+x)^ 

J (2^)2 q2(q2+^2) \ j 



Since Mi_o,4,C = (-Mi,o,4,c)*; we obtain 

f d^TJ f k d^o Uj^ 

Mi,o,3,c + Mi,oAC - -8LT gUjctatc ^L^. \JiP)? 



{2tt)^2E' " J (2^)32w (27r)2 q2(q2+^2-) 
fc2 / sm{^Ly 



APPENDIX F: COMPUTATION OF DIAGRAMS Mi,o,3,h AND Mi,o,4,i, 



We will now calculate the cut diagrams Afi^o.3,fl and Afi_o.4,L, shown in the Fig. (5). We start with AIi^.s^r: 



19 



1^1 ,0,3, R 



1^1,0,4,1 




x*-xJ<2L 



I 



FIG. 5: Feynman diagrams Mi,o,3,h and Mi,o,4,_l contributing to the radiative energy loss to first order in opacity, labeled in 
the same way as Fig. 1. 



1=0 

X A|^(a;3 - X2)v^{x3) A_+(x4 - X3) v^{x4) AZ^{x5 - 0:4) J(x5) 



where 



Here 



n 0(4+1 - 4)0(4 - 4) OC-L - (xi - xo)+) 



z^O 

00 /'OO 



dp^ (fpi f (Pk d^p d^q 4 



(27r)3 2p+7 (27r)32tj (27r)32£; (27r)4 

X (P2 +P3r^>A(g)b+Pl)^ (Pl +P2YPpa(k)(p+Pir I , 



» 5 3 

/ = / n^'^' n ^(4+1 - 4) - (xi - a:o)+) - 4) j(^o)j(x5) 

•' i=0 j=0 

X g-W3-(xi~Xo) p-iq-{x3-Xi) ^-ip2-{x2-Xi) ^-ipi-(x3-X2) ^-ip-(xi-X3) ^-ik-(xi~Xi) ^-ip4-{x5-X4,) 

= I J(p)|2 (27r)3 S{(p, -p- k)+)6\p; -p-k) (2^)3 5((p2 - Pi - k)+)S'(P2 - Pi - k) 
X (2^)3 6((pi + q - p4}S^Pi +q-p) (inf d{(pi -p- k)+)6^(Pi - p - k) h . 



(Fl) 



(F2) 



POO r^3 r^L, poo 

Jo Jo Jo Jo 

where we defined x'l = x\ — xq, x'2 = X2 — xi, x'^ = X3 ~ xi and x'^ — X4 — x^. By using 5 functions from Eqs. (F2) 
and (C5), we obtain in soft gluon soft rescattering approximation: 



(P3-p-k) = (Pi-p-k) ~(p2-k-pi) w 



(pi - pY 



q2 + 2kq 

P+ 



< C < |gUfe|,9^, leading to (pi -p) +q q- 



Ii then becomes 



= — (sinCL + ?:(l -cos^L)) / d4+(e-5(«"~«)<3'" -e-*?"^''). 
? Jo 



(F4) 
(F5) 

(F6) 
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Finally, by using Eqs. (F2) and (F6) Eq. (Fl) becomes 



-^•^1,0,3,7?, = 2i(7 tatctat 



[2'Kf2E 



\J{P)\ 



d^k 



fc2 



2Im 



1 1 



sinfL 1 — cos^L 

+ i 



oc 

X / dx'^{e' 



Note that A/i.o,4,L = (-^^i,o,3,ii)*, leading to 



(F7) 



-^^l,0,3,fl + ^^^1,0,4,L — ^Lg tatctcta 



d^p 



d^k 



fc2 



{27t)^2E' 7 (27r)32tj (A;2 

2 '1 1 , , 

1 - 



1 — cos 



where y = Lets define 



1 \ 



dy(sin(q" 



1-i 



Oy-sm{q-y))}, (F8) 



(F9) 



In high temperature limit F{qo, q^, q^) is even function of go a-nd q^. Then ^^^"g^^" '"^ 
tion of qz, leading to 



is an odd func- 



dS 



Fiqo,qi,q') (<5(go - Qz - - S{q2 - = ^ 



d^q dF{qo,ql,q^) 



(2^)4^^^"'"^'"'^^"" "V.. V (2^)3 dqo 

We will now compute the second part of the integral in Eq.(F8): 



<3o=<?2 



= 0. 



j^;^ F{qa,ql,q'^) J dyi^sm{qo ~ q^ - - sm{qQ - qz)y 



To do this we first concentrate on 

/•oo ^4 

''Jo (4 

OO /'OO 



''Jo (2.r 



F{qo,qz,q ) sin((qo - - Oy) 



^"(90, 9^, 9^) ( sm{qoy) cos{{qz - Ov) ~ cos{qoy) sin(((7^ - ^)y) 



cos((g^-Oy) / dqoF{qn,qz,q ) sin(goy) 



/ '^^ / '^{27t)^ cos(goy) y" dqzg{qo,ql,q^) svti{{qz-i)y)- 



(FIG) 



(Fll) 



Since ^ <C (?z, |<7| we can assume (/^ — ^ w q^. Therefore sm{{qz — ^)y) w sinq^y, leading to / dq^ F{qo, ql, q^) s\Tv{{qz — 
Oy) ~ / dqz F{qo,q'^,q^) sin{qzy) = 0. Similarly, for the second part of the integral in Eq. (FIG) we also obtain G, 
which finally leads to 



Ml,0,3,fl + Ml HAL « . 



(F12) 
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J(Xo) 




x|-xJ<2L 



FIG. 6: Feynman diagrams Mi^q.s.h and Mi,o,6,-L contributing to the radiative energy loss to first order in opacity, labeled in 
the same way as Fig. 1. 



APPENDIX G: COMPUTATION OF DIAGRAMS Mi,o,5,r AND Mi,o,6,i 
We will now calculate the cut diagrams A/i^o,5,i? and Mi_o,6,l, shown in the Fig. (6). We start with A/i_o,5,i?: 

X v+{x3) D''S_^{x4 - X3) A_+(a;4 - 0:3) (2:4) Al_(a;5 - X4) J(a;5) 

X e{x+ - x+) 0{x+ - x+) e{x+ - x+) e{x+ - x+) e{2L - (xi - .to)+) 



00 roc 



gXt^Uuj^ y_^n (2,)3 2p,+ 7 (27r)32a.(2.)32i? 7 (2.) 

X {(P2 + P3rD>{q)ip, +P2r} {{p + piYPp.{k){p+pir} i 



(Gl) 



where 



5 

I ^ W d'^x^9{x+ - x+) 9{x+ - x+) d{x+ - x+) d{x+ - x+)e{2L - {x+ - x+)) 

i=0 

X g-*P3-(2;i-a;o) g-iq-{x2-xi) ^-ip2-[x2-xi) ^-ipi-{x3-X2) ^-ikixi-xs) ^-ip-[xi-X3) ^-ipi-(x^-Xi) 

= I J(p)|2 (27r)3 5{{p3 -p- k)+)6^{p3 p-k) (2^)3 S{ip2 +q-p- k)+)5^{p2 +q-p-k) 
X (2^)3 6{ipi -p- k)+)5^ipi p~k) (27r)3 5{{p^ ^p- k)+)5\pi -p-k)h, 



J{xo) Jixs) 



(G2) 



and where 

p2L poo poo poo 

Ii= dx'+e-i^P^~P~''')~'='^^ dx^ei^P^-P-^^''^'^ I dx'^ e'-^^P^-P-'''^' ""'^ <i.T2+e-^(P^+9-P-'=)"^2+ . (G3) 
Jo Jo Jo Jo 

Here we defined x[ ~ xi — xq, x'^ = X2 — xi, x'^ = X3 — X2, x'^ = x^ — x^. 

By using S functions from Eqs. (G2) and (C5) we obtain p^ = pj = and: 

(P3 - p - k)- = {pi-p- k)- = {pi-p- k)- = -e (G4) 

Eq. (G3) then becomes 



16L^sinei_^cos^L-l 



^1 = -jfi—TT^) ^y^-Mp2-p-k+,r „ 

where y = x'^jl^ and in the last step we used (p2 — p — fc + g)^ = ^ S, ~ for small ^ 



(G5) 
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Finally, by using Eqs. (G2) and (G5), Eq. (Gl) becomes 

r (pri f (Pk k 



2 1 „2 



1 



X [-JJ^~^^J—) I dye - (G6) 

Since Mi o,6,l = (Mi_o,5,i?)*! it is straightforward to obtain (note that, similarly to Appendix D, part 
will vanish under integration over q): 



Mi,o.5,i? + A/i.o,6.L = -4.LT gUjatctc J (2^)3 2^ 1-^(^)1 J 



)^2uj (2^)2 q^(q^+fi^) 



APPENDIX H: COMPUTATION OF DIAGRAMS Mi,i,i,c AND Mi,i,2,c 

In appendices (H) - (K) we present in some detail the calculation of the diagrams where only one end of the 
exchanged gluon q is attached to the heavy quark, i.e. one end is attached to the radiated gluon k and consequently 
one 3-gluon vertex is involved in the process. In this appendix, we start with the calculation of the diagrams shown 
in Fig. 7. 

We will first calculate the cut diagram A/i,i.i,c: 



Ml, 



1,1, c 



/ Yl dxi J(a;o) A++(a;i - xq) v+ (xi) DXI^' {x2 - xi) w+ ^^^.^ (a:2)A_+(a;3 - a;i) 

X D'^'^{X3 - X2) V~{X3)D^^{X4, - X2)V~{X4) Al_(x4 - X3) Al_(x5 - X4,)Jix5) 

X e{xf - x^)0{x+ - xt)e{2L - {x2 - xo)^)0{xj - xi)e{xi - x+)e{2L - (.T4 - 2:5)+) 

dp+d^p, dk+d^k, f d^k d^p dS ifcba,,, 

Jo lX{2^f2pt{2nY2kt J (2^)3 2c. (27r)3 2£; (2^)4^ ^ '^'''^ 
X {p + piYP^^pAki) {gP'PHki+k)P- +gP^-P^{q-k)P' +gP'P-{~ki - q)P') 
X Pp3,(fc) D>{q)iP + P2r{P2+P:ir I. (HI) 
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where 

/ = 



5 

/ Jl dx^e{x+ - x+)e{x+ - x+)9{2L - {X2 - xo)+)9{x+ - x+)9{x+ - x+)9{2L - {x^ - 2:5)+) 

X g-ipi{xi~xo) ^—iki{x2~Xi) ^-ip{x3-Xi) ^-ik{x3-X2) ^-iq{x4-X2) p-ip2{'J:4.-X3) ^-ip3{x5-Xi) J(^xo) J{xz) 

\J{p)\''{2^f6{{pi -p- ki) + )5\pi -p- fci)(27r)3j((fci - fc - q) + )6^ki k q) 

X i27rfSiip + k- P2)+)S\p + k- p2){27TfS{{p + k + q- P2)+)S\p + k+q - P2) h , (H2) 



and where 



h = dx+ e"^('^i-'^-'?)-4+ / ' dx'^ e-W^-p-ki)-x\+ 



■i{P+k-p2)^ x'+ I ^^1+ ^- ^(p+k+q-p3y x'+ 



J, J. 3 



} 



'^{pi - h-p)- ^-^{Pi -p-k-q)- -|(fci -k-q)- 

1 ^-i{pi-p-k-q)- L-i _p_ fc_ g)- g-»(fei-fc-g)-i _ 1 



\{Pi-ki-p)- -^{pi- p-k-q)- \- e-*(Pi-J'-'=-9)"-^ - 1 -^{ki-k-qf 



(H3) 



By applying 6'^ functions from Eq. (H2), and by using p^ — p — wc obtain 



Pi = P3 



p+ {p + k + q)+ 



which leads to 



_ _ {k+q)^ + AP If+m^ {k+qf+ml _ _ q^ + AP 

^ ~ ^ ' ^ {k + q)+ ^^^bT^' ^^'^> 



iP+k-P2)- ^ (H5) 



(derived for x = ^^^^.^ ^ and .t <C 1) and 



(derived for xi = « « a;, and a; < 1). 

By using Eqs. (H5) and (H6) and assuming (as in previous sections) that finite size effects are negligible for collisional 
contribution, wc obtain 



^1 = n.2^..u,.^.^2 , .A ^---7^)2^L6{q,-q.), (H7) 



-16x'^E+^ / e-"^^ - 1 

(fc2+x)((fc+9)2+x)V^"~zCi 

where we used S{{p + k + q — pi)-) « (5(^1 — qz). 
Lets now compute 

(pi+P2)^(p+P2)''(P2+P3rPppi(fci)Pp3.(fc)^P2a(9) [gp'^Hk + k^r^+gp^p^iq-kr +g'''H-ki-qr') 
« { (pi + P2 Ppp, {k,)ppHk){p+p2r}{{k+k, YD>^ [q) {P2+P3r} 

- { - - 1) /(^°) ^^^^ S ^^-(^diiM - }■ 
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(for more details see Eqs. (C1-C4) in [8]). Also, note that 



Finally, by using Eqs. (H2), (H7), (H8) and (H9), Eq. (HI) reduces to 



(H9) 



Ml, 



1,1, c — 



(fc2 + x)((fc+q)2 + x) I 



f2E 
sin C_L 



\J{P)\- 



(Pk (fq 



+ i 



{27t)^2uj (27r)2 q2(q2 + ^2) 
1 — cos CL ~i 



(HIO) 



Since Mi.i,2,c is a complex conjugate of A/i^i.i^c", one finally obtains 



M 



1,1, 1,C + -^^^1,1,2,0 



k ■ (fc+g) 



(27r)32£; 

sin 
1 ^ 



(27r)32w (27r)2 q^{q^ + ix^) 



(Hll) 



APPENDIX I: COMPUTATION OF DIAGRAMS Mi,i,2,h + Mi,i,i,i 
In this appendix we will calculate cut diagrams A/i i i.^ and Afi.i_2.i?,, shown in the Fig. (8). We start with i\/i,i,2.i?,: 



x|-xJ<2L 



1,1, 1,L 




x\-xt,<2L 



FIG. 8: Feynman diagrams A/i,i,i,l and Mi,i,2,fl contributing to the radiative energy loss to first order in opacity, labeled in 
the same way as Fig. 1. 



5 

Mi,i,2,fl = / n '^^^•^(a^o)At+(a;i-xo)w+(xi)i^+7^(a;3-xi)i;+^^^3(x3)A++(x2 

•J n 



5 

, ^ . 

1=0 

X -D'^+(x3 - a:;2) {X2) D''_^"{xi - X2) v„{xi) A_+(a;4 - X2) wZ_(a-5 - a;4) ^(2:5) 
X e{x+ - Xfi") e{xi - xt) e{2L - {xi - a;o)+) e{xt - xt) e{xt ~ xt) 

-A- dp+d'^p, dk+d^h f d'^k d^p d^q 4 
,fj (27r)3 2p+ (27r)3 2fc+ i (2^)3 2c. (27r)3 2ii; (2^)4 ^ ^ '"'^''^ 

X {pi +P2rip+Pir{p+P3r D>p^{q)p,p,{ki)Pp,4k) 

X c/''i''=' (fc + ki)P^ + gP^P-' (-fc - qY^ + gP^P^ {q - k{)P'' /, (II) 
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where 



I = Jfl dx,9{xt - x+)9{x^ - xt)0{2L - {xi ~ xq)+)0{x1; - xt)0{xt - x^) 



^ ^-iP2(xi-Xo) ^-iq{x3-Xi) ^-ipi(x2-Xi) ^-iki(x3-X2) ^-ik{x4,-X3) ^-ip{xi~X2) ^-ip3ix5-X4) J(^xo) J{x^) 

= I J(p)|2 {2TTfS{{p2 -p- k)+)6\p2 -p- k){2TTf5{{pi -p- h)+)5^{pi -p-ki) 

X (2^)35((fci + q - k) + )S\k, +q- fc)(27r)35((p3 - p - k) + )d^p3 - p - k) h , (12) 

and where 

h = dx'+ e-5(P=-P-'=)~^'i^ / dx'+ e^^P'-P-'''^^^^ 

Jo Jo 

/>oo 

X / dx'+ e-5(fci+9-fc)"4+ / dx'^ ^-i(Pi-p-kl)-x'+ (j3) 

Jo Jo 
In the last equation we defined x'l — xi — xq, x'2 = X2 ~ xi, x'^ = X'i — xi and X4 = 2:4 — 2:5. 

By applying 5'^ functions from Eq. (12), and by using p^ — we obtain P2=P'3~ (p+k)+ ■> leading to (in soft 

gluon soft rescattering approximation) 

{p2-p-ky = (p3 -p- fc)- = 

(Pi-p-fci)- «-C, (14) 

where in the second relation we used xi = (^^^^'1^^+ ~ {p+k)+ ~ ^' 
Ii then becomes 

^ 8Ll-e^ /-""^ ,+ f _.(,-_j)4+_^_.(,-+C-CKn. (15) 
Furthermore, similarly as in [8], for highly energetic jets 

{{p+Pir Pup^ ih) {k)iP+P3r}{{Pi +P2rD,,Mik + fco"} 
. {_.Mi^}{.n..a-|,/(,4.i„(-^--^)}. 

By using Eqs. (12), (15), (16) and (H9), Eq. (II) finally reduces to 

4 r , , 1 r, ,1 / d^p 2 f d^k d'^q 



{2tt)^2E' '"^^ J (27r)32w {2tt) 
■^'(1 - ^J/Wo) ^ 2Im ' 



{k^+x){{k-qY + x) ' q2-'--q2 W-U^iq) q^-^riq) 

Since Mi 1 1 is a complex conjugate of A/i 1 2,^, it is straightforward to obtain [y = xt^ /2) 



/d^v f d 

(2^0fe j (2^ 



4 



f2uj (27r) 

■^'(1 - ^)/('?o) — 2Im ' 



(fc2 + ;^)((fc_q)2 + ;^) ^ cj2^^v-^^2 V9'-nL(<z) g2„nT(<z) 
X {^^(%o-g.-e)-%o-9. + C-0) 



1 — cos^L 



dy| sin(go - ~ Ov - sin(go - 9^ + C - 0^} ■ (K) 



By applying the same procedure as in Appendix (F), we obtain that Mi i ix + Mi_i^2,R. — 0. 



26 



APPENDIX J: COMPUTATION OF DIAGRAMS Mi,i,3,c AND Mi,i,4,c 



In this appendix wc will calculate the cut diagrams Mi^i^^^c and iV/i_i_4^c shown in the Fig. (9). We start with 
Mi,i,3,c: 



'1,1,3,C 




J(X5) J(Xo) 

0- 




P2 XufJ X2,V 



X2-Xq<2L 



1,1,4,C 



X4, o- p3 



-i 



I 



xJ-xJ<2L 



FIG. 9: Feynman diagrams Mi,i,s,c and Mi,i,4,c contributing to the radiative energy loss to first order in opacity, labeled in 
the same way as Fig. 1. 



Mi,i,3,c = / Y[dx,J{xo)A+^{xi~xo)v+{xi)DX'i'"{x2-xi)v+p^p^{x2)A^+{x3~xi) 

i=0 

X D'^^^{x3 ~ a;2)w^(x3)Al_(a;4 - X3)D'^^{x4 - X2) v~ {x4)AZ_{x5 ~ X4)J{x5) 
(27r)^ 2p+ 



-00 JQ 



X 6l(x+ - x+)9{x+ - x+)e{2L - {x2 - xq)+)9{x+ - x+)6'(x+ - x+)e{2L - (.T3 - X5)+) 



oc poc 



(27r)3 2p'+ (2^)3 2k+ J (27r)3 2c^ (27r)3 2E (2^)4 ^ '"^ 

X {gP'P'{ki+ky^ +gP'-P-'{q^kY' +gP'P'-{-ki-qY') I, (Jl) 



where 



5 

/ = / n - x+)9{x+ - x+)e{2L - {x2 - xo)+)e{x+ - x+)e{x+ - x+)e{2L ~ {x^ - X5)+) J{xo) J{x^) 

1=0 

■^g-iPl{^l--Co) ^-ikl{x2~Xi) ^~ip{x3-Xi) ^-iq{x3-X2) ^~ip2ix4-X3) ^-ip3{x5-X4) ^-ik{x4-X2) 

= I J(p)|2 {27tY5{{p, -p- ki)+)5\pi - p - k){2^f5{{ki - fc - g)+)52(fci - fe - 9) 

X {2T:f5{{p2 ~p- q)+)5\p2 -P- q)i2nfS{{p3 - P2 ^ k)+)5\p3 -p2-k)h, (J2) 

and where 

p2L r2x'^ 
h = dx'+ e-5(fei-fc-9)"4+ / dx'Y e-^^Pi-P-'^'i^^^'i"' 



2L r'^^'3^ 

dx'+ e-*(f^-P-«)"^3+ / dx'Y e" 5 (J3) 
Jo 

Here we defined x'^ = Xi — Xq, x'^ ~ X2 — xq, x'^ = x^ — x^ and x'^ = x^ — x^. 
By applying 5^ functions from Eq. (J2), and by using p^ = ^' '^^^ we obtain 

(fc+g)^ + AP _ _ + RP _ _ _ _ AP 

= K = ^-^P^^, (J4) 



k+ ' ^ (fc + g) 
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which leads to 

(pi - p - fci)" = -C, 
(derived for xi = j^^j^ ~ j;^^ « a: < 1) and 

-P2- k)~ = 

(derived for X2 = «^ ~ I^JTIF ~ ^ « !)■ 

By using Eqs. (J4)-(J6), Eq. (J3) becomes 



(J5) 



(J6) 



^ _ 16 r4sin^(pi -p- fc- g)-L/2 / {pi-p-k-q)' e-'^^^i i_i 
^ ~ C^l ((pi -p-fc-g)")^ V e-i(pi~p-fc-9)-i _ 1 (ki-k-q)- 



4 sin^ (P2 - P - g) -^ /2^ {P2-p-q) 
((P2 - P-?)")^ 

(p2 - p-qy 



(fci - fc - 9)" 

g-i(fei-fc-ij)"L _ 2^ 



g-i(P2-p-9)-L _ 1 (ki-k-q)- 



g-i(P2-p-«)-i _ 1 pi-p-k-q)- 



)} 



27rL 



16 



^27rL <S(g- 



-iQL 



) -<5(pi-p-g)-)( 



-^(C-«)i 



-z(C-Oi 



0} 



1 - 



-i(C~Oi 



(J7) 



where, as in the previous sections, we assumed that finite size effects are negligible for collisional contribution, and ii 
the last step we used Eq. (A9). 

Analogously to calculations performed in [8] (see Appendix C in [8]) one has 

{p + Pir{p + P2nP2+P:ir P^pAh) O^^iq) Pp,,ik){gP^PHk + k^r + gP^P^iq - kr + g^^^'^i^k, ~ qr) 

« {p+PiT ^MPi(fci) ^r(fc) {P2+P3r[{k + k,r D>jq)(p+P2r] 

By using Eqs. (J2), (J7), (J8) and (H9), Eq. (Jl) finally reduces to 

d^p 



Ml 



1,3, c 



-2LTg*[ta,t,] [t,,ta] 
k ■ (k+q) 



{2ttY2E 



\J{pt 



(Pk (Pq /i^ 



(J8) 



(27r)32tj (27r)2 ^2(^2 ^ ^2) 



-iQL 



(fc2 + x)((fc+g)2 + x 

Since Mi.i^^c is a complex conjugate of Afi^i.a^c*, one finally obtains 

cPp 



(J9) 



fc • (k+q) 



(Pk (fc 



'{k^ + x){{k+qr + x) 



■)32w (27r)2 q^{q^+pfl) 
sin^i sin^L sin(C — ^)L' 



(JIO) 



APPENDIX K: COMPUTATION OF DIAGRAMS Mi,i,3,h , Mi,i,3,i , Mi,i,4,i AND Mi,i,4,h 



In this appendix we will calculate the cut diagrams Mi^i^^^n, Mi^i^^^l, Mi^i^^^n and Mi_i_4_i shown in the Fig. (10). 
We start with Mi 1 3 
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FIG. 10: Feynman diagrams Mi,i_3,_r, Mi,i,3_l, ^1,1,4,^1 and Mi^ia,l contributing to the radiative energy loss to first order in 
opacity, labeled in the same way as Fig. 1. 



5 

Mi,i^3,R = / Y[dx,J{xo)AX^{xi~XQ)v+{xi)DX^l'"{x2-xi)v+p.^p^{x2)AX+{x3~xi) 



i=0 

-+ 



X D\''_^'-' {x-i - xi)v'l{xz) D''J'^{x4 - X2)v„ {xi) A_+(a;4 ~ xs) A (2:5 - Xi)J{x^) 

X e{x+ - x+)o{x+ - x+)e{2L - {x2 - x^)+)e{x+ - x+)e{x+ - x+) 

^ loo io iX (2^)3 2p+ (27r)3 2fc+ j {2^Y 2c. (2^)3 2E {2^Y ^ ^ 

X {(pi + Ppp, (fci) (fc)(p + Ps)'^} {(fc + k,Y^ D>,iq){p + prr]l , (Kl) 
where we used that in soft gluon soft rescattering approximation (see Eqs. (C1-C4) in [8]) 

{Pi + P2r {p + PiYiP + Ps)'^ P^.p^ (kl) D>_M PpsAk) {g'^^Hki + kr- + gP-P-iq - kf^ + gP^P- {-k, - qf') 

^[{Pi+ P2T P^P^ (kl) PS' {k){p + PzY }{{k + k^r D>MiP + Pir} • (K2) 
In Eq. (Kl) / corresponds to 



/ = 



5 

/ n dx,9{xi - x+)9{x+ - x+)ei2L - (.T2 - xo)+)0(.T+ - x+)eix+ - xt) 

^ ,: — n 



X Q-iP2ixi-xo) ^-iki(x2-xi) ^-ipi(x3-xi) ^-iq{x3-X2) ^-ik{x4-X2) ^~ipix4-X3) ^-ip3{xB-X4^ J{xo) J{x^) 

I J(p)P {2tt)H{{p2 -pi- ki)+)6Hp2 -Pi- fci)(27r)35((fci - k + pi - p)+)5\ki - k + pi - p) 

X (2^)3<5((pi + q - p)+)5\pi +q p){2^f5{{p3 -p- k)+)S\p3 -p-k)h, (K3) 



where 



2L 



h = I dx'+ e-i^'''-''-P'-P^ "^S^+y dx'+ e~i^P^~P'-''''> <^ 

dx'+ e-^^Pi+'-P)"^'^"" / dx'^ ei(p^-P-kr A^ (K4) 
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By applying 5^ functions from Eq. (K3), and by using = ^' '^^^ we obtain in soft gluon soft rescattering 
approximation (see also (A9)) 



P2 =P3 



E+ ' 



{p2 -pi- = -c, 

{p3-p- ky = (p2-p- ky ^ , 

(fci - fc + pi -p)" = c - 

{q + Pi -pT ~qo- qz- 



By using Eq. (K5), Eq. (K4) becomes (y = ^) 



h 



-iq y 



dy. 



Similarly to previous appendices, Eq. (K2) is equal to 

{ (pi + Ppp, (fci ) (fc) (p + P2 r } { (fc + fci r , (g) + } 

k-{k+q) 



-4- 



1 - 



2Im 



Finally, by using Eqs. (K3), (K6)-(H9), Eq. (Kl) reduces to 



k ■ (k+q) 



d^p 



q'--IlL{q) q^-TlT{q) 
d^k d^q 



{2TTf2E 



\J{P)\ 



((fc+q)2+x)((fe2 + x) 



X e{ 1 



(27r)32cj (27r)4 



1 



Since Mi i 4 is a complex conjugate of Mi i^^ fj, one finally obtains 

d^p 



q'--nL{q) q^-Uriq) 

-'"'ydy. 



k ■ {k+q) /sin^L sin(C — C)-^ 

{{k+qf+x){k^ + x)^~lL {C-OL 

Similarly A/i_i_3_l + Afi,i_4_ij = Mi,i,3,fl + Mi^i^i^L leading to 



\Jip)\' 



d^k d'^q 



(27r)32tj (27r)2 q^{q^ + ^i"^) 



k ■ (k+q) 



d^k d^q 



(K5) 



(K6) 



(K7) 



(K8) 



(K9) 



(27r)32u; (27r)2 



sin^L sin(C — 
q2(q2 + ^2) ((fe+q)2 + ^)(fc2 + ^) V^Z (C-Oi 



) 



(KIO) 



APPENDIX L: COMPUTATION OF DIAGRAM Mi,2,c 

In the Appendices (L) and (M) we present in some detail the calculation of the diagrams where both ends of the 
exchanged gluon q are attached to the radiated gluon fc, i.e. there are two 3-gluon vertices involved in the process. 
In this appendix, we start with the calculation of the diagram Mi^2,c shown in Fig. 11. 
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FIG. 11: Feynman diagram Afi,2 contributing to the radiative energy loss to first order in opacity, labeled in the same way as 
Fig. 1. 



5 

Afi,2.c = / ]ldx,Jix„)A++{xi-xo)v+ix,)D+'l'''{x2-xi)vXp^_p.^ix2)DPj^'{x3^X2) 

i=0 

X D'^^^{x3 - X2) A_+(a;4 - a;i)w^^ ^^^3(0:3) DZ'L"'{x4 - xg) v~ (x^) Al_(x5 - X4)J{x5) 
X 9{x+ - x+)9{x+ - x+)9{2L - (xa - xo)+)0{x+ - x+)9{xt - x+)e{2L - {x^ - xs)^) 

- ^--^ dp+d'p, dk+d^h f d^k d^p A 4 , 
, fjf (27r)3 2p+ (27r)3 2fc+ J (27r)3 2c^ (27r)3 2^; (2^)4 ^ 

X {^{k + k,r- D>_,^{k + k^r) (^{p+p^rp^pAki)p^'''Hk)PpAk2){p+P2r)i, (li) 

where in the second step we only keep the dominant parts from triple gluon vertices (sec [8]). 
In Eq. (LI) / corresponds to 

5 

I = Yi d■x^e{x+ - x+)e{x+ - x+)9{2L - {x2 - xo)-^)0{2L - (xg - X5)+)6l(x^ - x|)6l(x+ - xt) 

X Q-Vl{^l-^o) Q-ikl{x2-Xi) ^-ik{x3-X2) ^~iq{x3-X2} ^-ip{x4,-Xi) ^~ik2ix4-X3) ^-ip2{x5-X4) J(^Xi)) J(X5) 

= I J(p)|2 (2^)3 5{{pi -p- k,)+)6^{p, p fci)(27r)35((fci -k- q)+)S^{k, - k - q) 

X (27r)35((fc2 - q - k)+)5\k2 q k){2^f5{{p2 -p- k2)+)d\p2 - p - fca) /i , (L2) 

where 



2L 



h = I dx'+e-i^'''-''-'^^ ""'^^ j dx'+ e-s^Pi-P-'^i) =^'1^ 
dx'+ ei'^^^-'i-''^' ""'^^ / dx't ei(.P^-p-k2)-x'+ 







Here x'^ = x^ ~ x^^ x'^ = xj — x,^, x^ = X;^ — x^ and x'^ = x\ — x^. Furthermore, from b functions in Eq. (L2) 
it follows 

P2=Pi, k2 = ki, {pi - ki - p)~ = {p2 - k2 - p)~ = -C ■ (L4) 
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After using Eq. (L4), Eq.(L3) becomes 



16 [Asin^ 



{{pi-p-k~q) L/2)-l ^ _ jpi-p-k-q) 
{{Pi - P - k - p)-y I e*(Pi-P-'=-9)"-^ - 1 



g-i(fci-fc-ij) L _ ^ 



ki — k — q 



(pi-p-k-q)- e'(fei-fe-9) i _ 1-, Asm^{{ki-k-q)-L/2)-l 



g-j(pi-p-fc-9)"i - 1 ki-k-q 



} 



{{k,-k~q)-r 



16 r 

2ttL 



27TL5{ipi-p-k~qy){ 



1 - -— T + 27rL5{{ki -k-q)~) 



iCL 



-iCL 



16 r 



GAttL 



sin 
sin^L^ 



(l - 2 - p - fc - q)-) + <5((fci - fc - q)-) 



(L5) 



As in the previous sections, we here assumed that finite size efi'ects are neghgible for colhsional energy foss. Further- 
more, in the last step we used Eq. (A9), i.e. {pi — p — k — q)^ sa (fci — k — q)^ ~ qo ^ Qz- 



Similarly as in [8], for highly energetic jets 

(p + p^r Pt.p. (ki) P^'"' (fc) PaAk2){p + PiY) - (p + PiT P,.p. (fci) P"^"^ (k) p,Aki){p + PiY) 



{k + k,r,D>,^iq){k + k2r « fc+fc+^(l-^)/((?o)^2Im 



1 



1 



q^^nUq) q^-TlT{q) 



(L6) 



Finally, by using Eqs. (L2), (L5) and (L6), and after performing the same procedure as in the previous appendices, 
Eq. (LI) reduces to 



1,2.C 



dPp 



{2ttY2E 



\J{P)\' 



(Pk cPq 



P 



{k+qf 



{{k+qf+xf 



sin(CL) ■ 



(L7) 



APPENDIX M: COMPUTATION OF DIAGRAMS Mi,2,ii AND Mi,2,i 



In this appendix, we will calculate the cut diagrams Mi^2.R and Mi.2,l, which are shown in Fig. 12. We start with 




FIG. 12: Feynman diagram Afi,2 contributing to the radiative energy loss to first order in opacity, labeled in the same way as 
Fig. 1 
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5 

X dX+'^''{x3 -X2) A_+(x4 - xi)w+^^^3(.T3)i:)°;i+(x4 -X3)U~(X4) Al_(a;5 -^4) J(X5) 
x6l(a;+ - .T+)6l(a;+ - x+)e{x+ - x+)e{x+ ~ x+)9{2L - {x2 - .To)) 



00 nOO 



n 



(fik 



abdx 



J -00 t\ (27r)3 2pt (27r)3 2fc+ 7 (27r)3 2uj (27r)3 2i? (2^)4 

X [{ki + k2rD>,{q){k + k^r}[{p+p,rp^p,{k2)p^^p^{ki)Pp,4k){p+p2r}i, 

where in the second step we only keep the dominant parts from triple gluon vertices (see [8]). 
In Eq. (Ml) / corresponds to 

5 

I = Yi dx,0{x+ - x^)0{x+ - xt)e{2L - {X2 - xo)^)e{x^ - xt)e{xt - xj) 



I J(p)|2 (27r)3 5((pi ~p- k2) + )S^{p, ~p~ k2){2^fS{{k2 ~ k) + )S\k2 - k) 

X {2T:f5{{ki +q- k)+)5\ki +q- k){2T:f5{{p2 - k - p)+)5\p2 - p - k) h , 



J{xo)J{x5) 



(Ml) 



(M2) 



and where 



2L 

+ 



Jo 



dx'+ e-i^''^+'^-''^ / dx'4^ eiiP2-p-k) x'+ 



(M3) 



In the last equation, we defined x'^ = x^ — Xq , x'^ = — Xq, x'-^ — x^ — and x'^ ~ x'\ — x^. 

After applying 5 functions from Eq. (M2) and by using Eq. (C5), we obtain k~ = kl^ , which reduces the Eq. (M3) 



to 



2L 



dx'.^ 



'0 

16L 



e^^^i / e-^«^4 / dx'+ e~^'J ^ 



sin^L 1 — cos^L 
1 — h i- 



dye 



-iq y 



(M4) 



In the above equation, we also used (pi ^ k — p) = {p2 — p — k) ~ — ^, (fci + q — fc) w 5 and y ~ ^ 
Similarly as in [8], for highly energetic jets 



{p+pirp^,p,{k2)pp'''Hk)PaAk)ip+P2r) = -— 



[k + kiYD>{q){k + kiY « k+kte{l - ^)fiqo) ^ 2Im 



1 



q^-IlL{q) q^-nriq) 



(M5) 



Finally, by using Eqs. (M2), (M4) and (M5), and after performing the same procedure as in the previous appendices, 
Eq. (Ml) reduces to 



Mi^2,R = -4LTg''[t,,ta][ta,tc] J 



d^p 



{2TTf2E 



\J{P)\' 



d^k d^q 



fc2 



{k^+xr 



1 



sin^L . 1 — cos^L 

~1l 



+ 1 



X 9{l-^)-^2lm 



(27r)3 2w (27r)4 
dye-"^'y 

1 



q^^ULiq) q^~TlT{q) 



(M6) 
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Finally, by using that Mi.2,l is a complex conjugate of Mi,2,_r, and by performing the same procedure as in 
Eqs. (D24)-(D28), one obtains' 



/d V f d 



d^k d'^q 



)3 2w (27r)2 q2(g2 +^2) 



fc2 



sin^L 



(M7) 



APPENDIX N: COMPUTATION OF THE TADPOLE DIAGRAMS 

1. In this appendix wc calculate tadpole diagrams shown in Fig. (13), and show that they present negligible 
contribution to the radiative energy loss. 




P I P 

I I 

FIG. 13: Tadpole Feynman diagrams Mt,R and Mt,L, labeled in the same way as Fig. 1. 



r ^ 

Mt,R = ]ldx,J{xo)A+4x,-xo)v+{xi)DX'lP{x2-Xi)v+^^^{x2)D+'l_-{x2)D'L\^^^ 

i=0 

X A_+(2;3-.Ti)Al_(.T4-X3) J(X4) 9{x+ -x+)9{x+ ^x+)e{2L~{x2-XQ)+)9{x+ -x+), (Nl) 
where Wp^Arl^^z) = -ff^ fceafcebC^gpagxr - QpXQar - 9pT9a\) , leading to 



dxi J{xo) 

i=0 

X (-1) 



ipi(xi-xo) 



d^p 



{2TTf2E 



OC /"OO 



(2^)3 2p+^ 



(2^)3 2k-l 

-^D^^>(a)(-1) ( .-^Hx,-.,)pau,,. f dp+d^p2 

(2,)4^ (?)( l)y (2^)32^^ P (2^)3 2p+ 



e-'P-(-^-^-){ig{p + P2)M) 



X e{x+ - x+) 0{x+ - x+) 0{2L - {x2 - xo)+) eixj - x|) J(X4) 



d^fc d^p dS A g^P^^V^ rffc+rf2fc^ 



(27r)3 2tj (27r)3 2E (27r)4 11 (27r)3 2p+ (27r)3 2fc,+ 
X (p + Pi)pP''''{ki){2gp,g^r - 9px9ar - 9pr9ax) D-^>{q)P-''{k){p+p2), /, 



(N2) 
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where 

f ^ 

1=0 

X e{x+ - x^)e{x+ - x+)0{2L - {x2 - xi)+)e{x+ ~ xt) 

= I J(p)|2 (2^)3 6{{p, -p- h)+)6\p, p k,){2nf6{{h ~ k)+)S^{ki - k) 

X {2^f5{{p2-p~k)+)5\p2-p-k)h, (N3) 

and where 



h - 



2L 



Here x'^ — x^ — Xq , x'^ = — Xq and .i 



(N4) 



x^ . Also, by using S functions from Eq. (N3) and Eq. (C5) it 



foUows 

P2=pi] k = ki; {pi - p - k)~ ^ {p2 ~ p - k)' 
After using relations from (N5), Eq.(N4 ) becomes 



2L 



dx:+ 



sin^L 1 — cos^L 

1 — h «■ 



(N5) 



(N6) 



By using Eqs. (N6) and (N3), Eq. (N2) becomes 



MtR = 8g^L[ta,t,][t„ta] 



d^p 



d^k d-^Q 1 1 



(27r)3 2uj (2^)4 {2E+f 2fc+ 



X {8p'^P^p{k)PS{k)p''D>\q) - 8p^^Ppp{k)DP'^>{q)P,,{k)p'') 



x^E+^ 
^ (fc2 + 

To proceed further, lets first calculate 



sin£L 1 — cos£L 
1 — h i 



(N7) 



We here use 



leading to 



fc2 

x"^ 



pf^P^^ik)PP{k)p'' = ---^; P^{q)^2; Q^^{q) = 1 



D>\q)^e{l 
Furthermore, in Coulomb gauge 

leading to 



(l + /(g°))2Im 



1 _ 9(1 

q2 



q^-Uriq) q^-ULiq) 



Ppp{k)QP%q)P„x{k) = G., 



(N9) 



(NIO) 



(Nil) 



p^P^,p{k)DP''>{q)P„,{k)p'' = p^Ppp{k)PP%q)P,xp'' 



p2 I jpk)^ I (Pl)^ ^{pk){pq){qk) ^ {pk){qkf 



fc2 q^ 



{k^y^q^ 



1 9 2 1 J_ (^^) 



1 _ 9cL 

q2 



92_nT(q) g2_ni(g) 



(N12) 
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By using Eqs. (N8)-(N12), Eq.(N7) becomes 



\J{pt 



sin 1 — cos 

1 — h i 



{2tt)^2E' " J (27r)3 2w (27r)4 (fc2 + 
:+ J (2^)4 V go 



-Uriq) 



Since Mtx is a complex conjugate of Mt^n, one finally obtains 



Mt,R + Mt,L = Sir/ [ta, t,][tc, ta] 



(27r)3 2E 



fc2 



(27r)3 2uj (fc2 + Af.T2 + 77152)2 



1 - 



sin 



where 



/* = 



d^q 

to J (27r)4 



1 - 



(qkf 



nriq) 



(N13) 



(N14) 



(N15) 



Finally, we want to show that It <C J ^^v{q'^) (where viq^) = ^r^^q:^^), leading to the conclusion that tadpole 
contribution is negligible to the first order in opacity radiative energy loss. To do this, we first observe that < 
< 1, leading to 



where 



d*q 

to J (27r)4 
where we defined ?/ = and 



0<lf < It. 



l-4)-2Im( 2 n r ^ 
q2/ go - ^T{q) 



7 rf^2^ 

47r UJ 



T /-2n _ 4 /"^ / 1 \ 

Jo X \q2(l - x2) + nT(x)/ 



(N16) 



(N17) 



(N18) 
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FIG. 14; Ratio ' is shown as a function of momentum |q|, for Debye mass ^ = 0.5. 

Furthermore, we want to show that It ^ J ^^^"'^(q^)- To do this, we have to prove that Jt,max{^^) is comparable 
with w(q^). We numerically confirmed that ^ 1; as demonstrated in Fig. 14 for typical value of Debye 

mass jj. = 0.5. That is, in Fig. 14 we see that the absolute values of the ratios are notably smaller than 1. We also 



Iql [GeV] 
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checked that the same conclusion is vaUd independently on the value of Debye mass. Having in mind that ^ 1, 

and by using It ^ J ^"^^(Q^); it becomes evident that h ^ j ^%fv{c^), which leads to the conclusion that tadpoles 
present a negligible contribution to the 1"* order in opacity radiative energy loss. 
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